par(mfrow=c(1,2))
seq=seq(0.005,0.995,0.01)

K1=NX1+2.576*NY1
L1=NX1-2.576*NY1;L1[L1<0]=0
K2=NX2+2.576*NY2
L2=NX2-2.576*NY2;L2[L2<0]=0

stap=0.01;seq=seq(.005,.995,.01)
D1=dnorm(seq,0.40,0.04);D1=D1/(sum(D1)*stap)
D2=dnorm(seq,0.70,0.10);D2=D2/(sum(D2)*stap)
G1=(D1+D2)/2
G2=dnorm(seq,0.45,0.15);G2=G2/(sum(G2)*stap)
plot(seq,NX1,col='magenta',type='l',lwd=2,ylim=c(0,8),axes=FALSE,
main=c(paste('6A. Respondent Profiles: Reference Cliffs Delta= .29'),
paste('Estimated Profile-Based Cliffs Delta = .16'),
paste('Its 99% Confidence Interval = (.10, .21)'),
paste('Group 1: 5 Items, Group 2: 5 Items')),
xlab='Probability of Behaviour',ylab='Plausibility')
axis(side=1,at=seq(0,1,0.10));axis(side=2,at=seq(0,8,2),las=1)
lines(seq,NX2,col='green3',lwd=2)
lines(seq,G1,col='blue',lwd=2,lty=2)
lines(seq,G2,col='red',lwd=2,lty=2)
lines(seq,K1,col='magenta',lty=2)
lines(seq,L1,col='magenta',lty=2)
lines(seq,K2,col='green3',lty=2)
lines(seq,L2,col='green3',lty=2)
legend('topright',bty='n',cex=.9,
c('Estimated Profile Group 1','Reference Profile Group 1',
'Estimated Profile Group 2','Reference Profile Group 2'),
col=c('magenta','blue','green3','red'),lty=c(1,2,1,2))

K1=NX3+2.576*NY3
L1=NX3-2.576*NY3;L1[L1<0]=0
K2=NX4+2.576*NY4
L2=NX4-2.576*NY4;L2[L2<0]=0

stap=0.01;seq=seq(.005,.995,.01)
D1=dnorm(seq,0.40,0.04);D1=D1/(sum(D1)*stap)
D2=dnorm(seq,0.70,0.10);D2=D2/(sum(D2)*stap)
G1=(D1+D2)/2
G2=dnorm(seq,0.45,0.15);G2=G2/(sum(G2)*stap)
plot(seq,NX3,col='magenta',type='l',lwd=2,ylim=c(0,8),axes=FALSE,
main=c(paste('6B. Respondent Profiles: Reference Cliffs Delta= .29'),
paste('Estimated Profile-Based Cliffs Delta = .24'),
paste('Its 99% Confidence Interval = (.19, .28)'),
paste('Group 1: 20 Items, Group 2: 20 Items')),
xlab='Probability of Behaviour',ylab='Plausibility')
axis(side=1,at=seq(0,1,0.10));axis(side=2,at=seq(0,8,2),las=1)
lines(seq,NX4,col='green3',lwd=2)
lines(seq,G1,col='blue',lwd=2,lty=2)
lines(seq,G2,col='red',lwd=2,lty=2)
lines(seq,K1,col='magenta',lty=2)
lines(seq,L1,col='magenta',lty=2)
lines(seq,K2,col='green3',lty=2)
lines(seq,L2,col='green3',lty=2)
legend('topright',bty='n',cex=.9,
c('Estimated Profile Group 1','Reference Profile Group 1',
'Estimated Profile Group 2','Reference Profile Group 2'),
col=c('magenta','blue','green3','red'),lty=c(1,2,1,2))


par(mfrow=c(1,2))

K1=NX5+2.576*NY5
L1=NX5-2.576*NY5;L1[L1<0]=0
K2=NX6+2.576*NY6
L2=NX6-2.576*NY6;L2[L2<0]=0

stap=0.01;seq=seq(.005,.995,.01)
D1=dnorm(seq,0.40,0.04);D1=D1/(sum(D1)*stap)
D2=dnorm(seq,0.70,0.10);D2=D2/(sum(D2)*stap)
G1=(D1+D2)/2
G2=dnorm(seq,0.45,0.15);G2=G2/(sum(G2)*stap)
plot(seq,NX5,col='magenta',type='l',lwd=2,ylim=c(0,8),axes=FALSE,
main=c(paste('6C. Respondent Profiles: Reference Cliffs Delta= .29'),
paste('Estimated Profile-Based Cliffs Delta = .27'),
paste('Its 99% Confidence Interval = (.23, .30)'),
paste('Group 1: 50 Items, Group 2: 50 Items')),
xlab='Probability of Behaviour',ylab='Plausibility')
axis(side=1,at=seq(0,1,0.10));axis(side=2,at=seq(0,8,2),las=1)
lines(seq,NX6,col='green3',lwd=2)
lines(seq,G1,col='blue',lwd=2,lty=2)
lines(seq,G2,col='red',lwd=2,lty=2)
lines(seq,K1,col='magenta',lty=2)
lines(seq,L1,col='magenta',lty=2)
lines(seq,K2,col='green3',lty=2)
lines(seq,L2,col='green3',lty=2)
legend('topright',bty='n',cex=.9,
c('Estimated Profile Group 1','Reference Profile Group 1',
'Estimated Profile Group 2','Reference Profile Group 2'),
col=c('magenta','blue','green3','red'),lty=c(1,2,1,2))

K1=NX7+2.576*NY7
L1=NX7-2.576*NY7;L1[L1<0]=0
K2=NX8+2.576*NY8
L2=NX8-2.576*NY8;L2[L2<0]=0

stap=0.01;seq=seq(.005,.995,.01)
D1=dnorm(seq,0.40,0.04);D1=D1/(sum(D1)*stap)
D2=dnorm(seq,0.70,0.10);D2=D2/(sum(D2)*stap)
G1=(D1+D2)/2
G2=dnorm(seq,0.45,0.15);G2=G2/(sum(G2)*stap)
plot(seq,NX7,col='magenta',type='l',lwd=2,ylim=c(0,8),axes=FALSE,
main=c(paste('6D. Respondent Profiles: Reference Cliffs Delta= .29'),
paste('Estimated Profile-Based Cliffs Delta = .28'),
paste('Its 99% Confidence Interval = (.26, .30)'),
paste('Group 1: 200 Items, Group 2: 200 Items')),
xlab='Probability of Behaviour',ylab='Plausibility')
axis(side=1,at=seq(0,1,0.10));axis(side=2,at=seq(0,8,2),las=1)
lines(seq,NX8,col='green3',lwd=2)
lines(seq,G1,col='blue',lwd=2,lty=2)
lines(seq,G2,col='red',lwd=2,lty=2)
lines(seq,K1,col='magenta',lty=2)
lines(seq,L1,col='magenta',lty=2)
lines(seq,K2,col='green3',lty=2)
lines(seq,L2,col='green3',lty=2)
legend('topright',bty='n',cex=.9,
c('Estimated Profile Group 1','Reference Profile Group 1',
'Estimated Profile Group 2','Reference Profile Group 2'),
col=c('magenta','blue','green3','red'),lty=c(1,2,1,2))


